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Our current understanding of strongly correlated electron systems is based on a homogeneous 
framework. Here we take a step going beyond this paradigm by incorporating inhomogeneity from 
the beginning. Specifying to systems near the Mott metal-insulator transition, we propose a real 
space picture of itinerant electrons functioning in the fluctuating geometries bounded by interfaces 
between metallic and insulating regions. In 2+1-dimensions, the interfaces are closed bosonic strings, 
and we have a system of strings coupled to itinerant electrons. When the interface tension vanishes, 
the geometric fluctuations become critical, which gives rise to non-Fermi liquid behavior for the 
itinerant electrons. In particular, the poles of the fermion Green’s function can be converted to 
zeros, indicating the absence of propagating quasiparticles. Furthermore, the quantum geometric 
fluctuations mediate Cooper pairing among the itinerant electrons, indicating the intrinsic instability 
of electronic systems near the Mott transition. 


I. INTRODUCTION 

The emerging picture of strongly correlated electron 
systems is that they possess a multiplicity of nearly de¬ 
generate ground states m- These systems are quali¬ 
tatively different from those with a single stable ground 
state, like Landau Fermi liquids and symmetry broken 
phases. A small change in the external parameters, 
e.g. pressure, doping, magnetic field, substantially al¬ 
ters the macroscopic properties of these systems. A 
grand challenge of modern condensed matter physics is 
to build a theoretical framework for such systems with 
large ground state degeneracy. When the excitation en¬ 
ergy is higher than the energy barrier between different 
ground states, the system is intrinsically inhomogeneous, 
with different ground states coexisting and fluctuating in 
real space [HUS]. The current paradigm of many body 
systems based on a homogeneous framework, with in¬ 
homogeneities treated as perturbations, is not adequate 
for the description of such intrinsically inhomogeneous 
systems. Here, we propose a framework that explicitly 
incorporates ground state degeneracy and inhomogene¬ 
ity from the begining. The basis of our formalism is 
the quantum geometric fluctuations associated with in¬ 
terfaces between different degenerate ground states. The 
study of emergent phenomena at interfaces has recently 
attracted much attention CHini, which provides further 
motivations for our work. 

For concreteness, we consider the Mott metal-insulator 
transition (MIT) |llll2j , which is prototypical of strongly 
correlated electron systems. At the MIT, the metallic 
phase and the insulating phase are degenerate in free en¬ 
ergy. In a coarse-grained picture, each coarse-grained 
region can be treated as in either the metallic or the in¬ 
sulating state (see Fig0. Let us focus on the charge 
degrees of freedom. Then the insulating regions can be 



FIG. I: (Golor online) Real space configuration of coexist¬ 
ing metallic (M) and insulating (I) regions, separated by the 
interfaces. 


treated as inert, e.g. not contributing to transport or 
Cooper pairing. The metallic regions are the active de¬ 
grees of freedom. A crucial observation is that the metal¬ 
lic regions are bounded by fluctuating boundaries, and 
these boundaries are themselves active degrees of free¬ 
dom. Then our basic picture of a Mott MIT consists 
of regions of itinerant electrons coupled to the fluctu¬ 
ating boundaries. Such a picture is fundamentally dif¬ 
ferent from usual condensed matter systems, which are 
generally defined on a rigid background. By contrast, 
in our formalism, the electronic system near the MIT 
constitutes a physical system operating on a soft back¬ 
ground with quantum mechanically fluctuating geometry 
and topology. Thermally fluctuating geometries, which 
we draw analogy with at various points, have been widely 
explored in soft condensed matter physics USHUj. How¬ 
ever, the quantum counterpart, to the best of our knowl¬ 
edge, was previously encounted only in quantum gravity 
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m- The closest example in condensed matter physics is 
probably the idea of fluctuating stripes |20] . 

Such quantum geometric fluctuations have far-reaching 
implications for the response of a system near a second- 
order or weakly first-order MIT, for which the geometric 
fluctuations become (nearly) critical. The relevant phys¬ 
ical systems include cuprates mm, organic conductors 
|22l23j , iron pnictides |24l25j and heavy fermions |26l27j . 
These systems display in their phase diagram a metal¬ 
lic region with anomalous single-particle and transport 
properties vastly different from that of Landau Fermi liq¬ 
uids. Such strange metal or non-Fermi liquid phases are 
usually associated with abrupt changes of the Fermi sur¬ 
face, which signals an itinerant to localized transition for 
certain degrees of freedom. Strange metal phases are be¬ 
lieved by many to be key to the understanding of high 
temperature superconductivity EIEHEg. Indeed we find 
that the geometric fluctuations of the interface lead to 
anomalous scaling in the fermion self-energy. Further¬ 
more, due to inversion symmetry breaking around the 
interface, couplings not allowed in a homogeneous en¬ 
vironment are enabled EOEI], which induce attractive 
interactions among the fermions and promote pairing. 


II. GLOBAL PHASE DIAGRAM 

A central object in our consideration is the interface, 
whose morphology gives rise to a new dimension in the 
global phase diagram. Let us start by considering a first 
order MIT (see Figj2}i). Close to the transition point, 
the system consists of mixtures of the metallic and in¬ 
sulating phases, separated by the interface. The length 
scales associated with the fluctuations of the interface are 
much larger than typical microscopic scales such as the 
interaction range and inter-carrier spacing. Such length 
scale separation makes it difficult to treat the interfa¬ 
cial degrees of freedom within a microscopic Hamilto¬ 
nian. Hence we will use a coarse-grained description, 
where one represents the two phases as domains of ferro- 
magnetically ordered Ising spins: spin up for the metallic 
phase, and spin down for the insulating phase [32]. Con¬ 
centrating on the behavior near the interface, one can 
further integrate out the bulk Ising degrees of freedom 
to obtain an effective Hamiltonian for the interface. In 
d-dimensions, the Hamiltonian is of the general form 

[IS] 


ns = j dS (a + + kK^ , ( 1 ) 

with the d — 1-dimensional area element d5', the interface 
tension cr, the mean curvature id, the Gaussian curvature 
dd, and the corresponding bending rigidities k and R. 

As one approaches the critical endpoint, the interface 
tension vanishes as cr(T) ^ (1 — T/Tc)^^ with the criti¬ 
cal exponent fi (Widom scaling). The zero temperature 
interface tension cr(T = 0) varies from system to system, 
forming a three dimensional phase diagram (FigJ^). As 



FIG. 2: (Color online) (a) Schematic phase diagram of a first 
order MIT. The horizontal axis denotes the external parame¬ 
ter that tunes the system through the transition, the vertical 
axis denotes the temperature, and the red dot is the criti¬ 
cal end point, (b) The three dimensional global phase dia¬ 
gram for the metal insulator transiton, including a new axis 
representing the interface tension a = a {T — Q). As a de¬ 
creases, the critical endpoint is driven to lower temperatures 
(red solid line), and becomes a QCP at a = 0 and T = 0. 
(c) The three dimensional phase diagram projected to the 
tension-temperature plane. The droplet phase with a > 0 
corresponds to the shaded region in (b). The sponge phase 
is controlled by the tensionless QCP. (d) Typical interface 
configuration in the sponge phase, with characteristic length 
scale 


cr(T = 0) decreases, the temperature Tc of the critical 
endpoint decreases. When cr(T = 0) vanishes, Tc also 
vanishes, and there is a continuous phase transiton at 
zero temperature, which forms a quantum critical point 
(QCP). Projecting the three dimensional phase diagram 
to the tension-temperature plane, one obtains the phase 
diagram shown in Fig. H- Such a phase diagram has 
been extensively studied in chemisty, biology and soft 
condensed matter physics in the context of surfactant sys¬ 
tems prHTS] . When the surface tension is positive and 
large, it is energetically favorable to minimize the sur¬ 
face area, and the system is in the droplet phase, with 
droplets of the minority phase residing in the majority 
phase. When the tension is negative and large, with a 
positive bending rigidity n stabilizing the system, it is 
energetically favorable to maximize the area of the inter¬ 
face, but at the same time minimize the interface curva¬ 
ture. The system then consists of alternating layers of 
metallic and insulating phases, and the system is in the 
lamellar (for d = 3) or stripe (for d = 2) phase. 

The droplet phase and the lamellar/stripe phase are 
separated by the QCP. In the quantum critical region, 
the system consists of bicontinuous percolating domains 
of the metallic and insulating regions (see FigJ^), and 
this phase is usually termed the sponge phase in the study 
of surfactant systems [IMIHI- Mapping the system to 
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Ising spins, the droplet phase corresponds to the ferro- 
magnetically ordered phase, and the sponge phase to the 
paramagnetic phase that restores the Z 2 symmetry. The 
sponge phase can also be reached by dynamically melting 
the stripes [20] , 


III. THE MODEL 


We consider systems near a second-order or weakly 
first-order MIT, which are then in the sponge phase at fi¬ 
nite temperatures. As a first step, we ignore the effect of 
disorder, which can locally pin the interface fluctuations. 
Since the fermions move much faster than the interface, 
there is a separation of time scales. With the character¬ 
istic length scale of the sponge structure, the char¬ 
acteristic velocity of fermions, the Fermi velocity the 
characteristic velocity of the interface, the sound veloc¬ 
ity 'T’s, and the condition vp ^ Vs^ there are two char¬ 
acteristic time scales: and At time scales 

t < assuming the fermion mean free path to be 

much larger than the fermion propagation is ballis¬ 
tic. At time scales ^k/vf < t < ^k/vs^ the interface can 
be treated as quasi-static, and the fermions propagate in 
a random porous media (see Fig. |2i). This corresponds 
to the region of Knudsen diffusion 33] , where the domi¬ 
nant scattering process is the scattering of fermions with 
the boundary. In this region, fermions are largely con¬ 
fined to local cages with typical size and transport 
occurs through tunneling between neighboring cages. At 
time scales t > ixhsi the interface fluctuates strongly, 
fermions experience the averaged effect of such fluctua¬ 
tions, and the system is on average homogeneous. 

To capture the long time scale physics, we proceed in 
the spirit of Hertz-Millis |34l35j to describe the system 
in terms of fermions coupled to the bosonic collective 
modes, here the quantum fluctuations of the interfacial 
geometry. The partition function of the system involves 
a summation over different fermion field configurations 
and a summation over different geometries: 


Z = 




( 2 ) 

where V represents the real space configurations of the 
metallic regions, and dV their boundaries, Ss the in¬ 
terface action, Se the fermion action, and 5int the cou¬ 
pling between fermions and the interface. We treat the 
fermions, without coupling to the interface, as a renor¬ 
malized Fermi gas with the action Se = j^dfirdtCe^ 
where 


Ce=^\r,t)(^^^ +^ + (3) 

with effective mass m*, and effective chemical potential 
/i*. Focusing here on the charge degrees of freedom, we 
have suppressed the spin index. 




FIG. 3: (Color online) Typical interface configurations con¬ 
tributing to the bosonic string correlator {0{Xa)0{Xb)) (a) 
and {0{Xb)0{Xc)) (b). 


A. Modeling the interface 


When the interface is classical, its action takes the form 
Ss = 1-islT. Near the MIT where quantum fluctuations 
are important, we also include the temporal fluctuations 
of the interface in the action. Let us now specify to d = 2, 
for the following reasons: (1) simplicity of theoretical 
treatment, (2) prominent MIT materials, e.g. cuprates, 
organic conductors, indeed have layered structure. In this 
case, the fluctuating interfaces are closed bosonic strings. 
We use with a = 1,2 to represent the worldsheet co¬ 
ordinates of the string, and with /i = 0,l,'' - ,D — 1 
its space-time coordinates. Here D = d + 1, and in our 
case d = 2. It is also useful to bear in mind that an Eu¬ 
clidean quantum field theory in d-dimensional space at 
T = 0 can be mapped to a (d+l)-dimensional classical 
field theory. Hence the quantum string action in (2+1)- 
dimensional space-time can be mapped to the classical 
membrane Hamiltonian in 3-dimensions (taking d = 3 
in Eq. 0 ). When the string tension a vanishes, the 
string worldsheets form sponge-like structures in 2+1- 
dimensional space-time. 

The leading order term in the string action is the ten- 
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sion term, named the Nambu-Goto action m 

=(7 J (f'iy/g, (4) 

with g = det^a 5 , and gab = daX^dbX^ is the induced 
metric. 

Since we consider the system near a QCP associated 
with vanishing tension, higher order terms in the string 
action need to be included. Hence the string theory 
for interface fluctuations near the QCP goes beyond the 
standard textbook string theory [19] which involves only 
the Nambu-Goto action. The higher order terms are as¬ 
sociated with the curvature of the worldsheet (see for 
example |13ll4j and references therein). The local cur¬ 
vature properties of a surface are encoded in the two 
principal curvatures i^i, K 2 , which are the inverse of the 
principal curvature radii. The product of the two princi¬ 
pal curvatures is the Gaussian curvature K = K 1 K 2 , and 
their average is the mean curvature H = {Ki + i^2)/2. 

The mean curvature term (Polyakov-Kleinert action) 
is of the form PHOT] 

f I (5) 

In terms of the string coordinates this term can be 
written as 

sp = (6) 

where AgX^ = -^dai-^g^-^dbXg). 

The third term associated with the Gaussian curvature 
K is the Einstein-Hilbert action of the worldsheet field 
theory 


) =kJ d^^^K, (7) 

since the Gaussian curvature K is related to the scalar 
curvature Rhy K = Rl2. This term is a total divergence. 
The Gauss-Bonnet theorem relates such a term to the 
topology of the surface 

I d^^^K = 27rx, (8) 

where y is the Euler characteristic of the surface. Thus 
this term is only relevant when the topology of the world- 
sheet changes. 

It is quite useful to write the action in the Monge rep¬ 
resentation, Xq = = Vgt^ Xi = (^2 = X 2 = 

where Vs is the sound velocity of the interface (see Eig. 

t . The area element is then dS = d^XyT^TXVhJ^, with 
X = dXodXi, and 


Vh = 


dh dh \ 


f 1 dh dh\ 
dt^ dx ) 


(9) 


The low energy effective action (usually named the cap¬ 
illary wave action, see e.g. [38] and references therein) is 


obtained by expanding the action in powers of field gra¬ 
dients, assuming that the variation of h in space-time is 
smooth. The Nambu-Goto action then becomes 


5?) 



l + \{Vhf 


( 10 ) 


In the Monge gauge, the mean curvature reads 


H = V ■ 


Xh 

yi + (vh)2' 


(11) 


Hence the Polyakov-Kleinert term can be expanded to 
give 

S(^) = ^jd^X{X^h)\ ( 12 ) 

The Gaussian curvature term is associated with the 
global topology change, and does not enter the capillary 
wave action. In summary, the capillary wave action reads 
= f dxdtVgCcw^ with the Lagrangian 


a I I'dh\^ 
= 2 



K / I d^h 
2 dt‘^ 


d^h \ 
dx‘^ ) 


2 

(13) 


Luttinger liquid perspective 

In the above consideration, the interface fluctuations 
are treated as noninteracting. However, since the inter¬ 
face fluctuations are confined to Id, the effect of inter¬ 
actions can be dramatic. In fact, in Id, both the inter¬ 
acting bosons and fermions are governed by Luttinger 
liquid physics at low energy [39]. We can use a general 
scaling form for the capillary wave correlator, assuming 
the interactions have driven the interface fluctuations to 
the Luttinger liquid fixed point. An equivalent way to 
take into account such Luttinger liquid physics is to map 
the quantum mechanically fluctuating interfaces in 2+1- 
d to random surfaces in 3d Euclidean space, and consider 
random Gaussian surfaces with a general probability dis¬ 
tribution where m 

Sew = ^ J dqduiiq^ (14) 

with the stiffness H, and the roughness exponent 0 < C ^ 
I. One geometric measure of such random surfaces is the 
fractal dimension Df of the contours with equal height 
h(x, t), which relates the length I of a contour to its radius 
r, with I . The fractal dimension 1 < Df < d = 2 

is related to the roughness exponent by [40] 

Df = (15) 

We note that fractal distributions of dopants have been 
observed in cuprates, and they are correlated with the 
enhancement of the superconducting transition temper¬ 
ature [41]. 
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B. Fermion-interface coupling 

We now specify the form of the coupling between the 
fermions and the interface. Similar to electron-phonon 
coupling in layered systems with broken inversion sym¬ 
metry, e.g. in the presence of buckling [42], we expect a 
local electric field arising from the asymmetric environ¬ 
ment around the interface, with its direction perpendic¬ 
ular to the interface ISUEI]. The interface thus possesses 
electric dipole charges, and the local segments of the 
strings can be thought of as oppositely charged parallel 
plates. The formation of such electric dipoles at the inter¬ 
face is widely observed in various heterostructures, e.g. 
metal-semiconductor junction, LaAlOs/SrTiOs [434146] . 
In the Monge gauge, for a static interface, the local elec¬ 
tric field gives rise to the coupling Si^t = J dxdtVsCint^ 
where 

£int = ( 16 ) 

Consider a covariant generalization of such a coupling 
by introducing the antisymmetric local electromagnetic 
field tensor Here we assume the magnitude of 

the electric field to be fixed. The fermion density p = 
then couples to the worldsheet operator 0[X] = 
F^^e^^daX^dbXjj. The interaction term thus reads 

>Si„t = xj ( 17 ) 

The combination plays the role of the Kalb- 

Ramond 2-form field familiar in string theory ng. 


IV. FERMION SELF-ENERGY FROM 
INTERFACE FLUCTUATIONS 

We now examine how the interface fluctuations influ¬ 
ence the dynamics of fermions. The leading order correc¬ 
tions to the fermion action are the self-energy effect and 
the density-denstiy interaction. The self-energy correc¬ 
tion is of the form 

55(1) ^ j J d^XB^ABi’HXAMXB), (18) 

with 

^AB = -X^{0{XA)0{XB))mXA)i’HXB)). (19) 

The induced interaction is of the form 

55 ( 2 ) = J d^XA J d^XBVABp{XA)p{XB), ( 20 ) 

with the fermion density p = and 

Vab = X^{0(Xa)0(Xb)). ( 21 ) 


So the effect of the interface fluctuations is encoded in 
the bosonic string correlator 

Vab = (O(Xa)O(Xb)} = [ VXO(XA)0(XB)e-^‘, 

( 22 ) 

which involves a summation over worldsheets that pass 
through Xab^ with the vertex operator O inserted at 

The full string theoretical calculation for such a rigid 
string theory is beyond the scope of the present paper, 
and we will proceed using capillary wave theory, which 
captures the low energy fluctuations of the interface. 
To calculate Vabi we employ the Monge gauge, where 
-^0 = = Vgt, Xi = ^2 = X 2 = h{t,x). Since 

on average the system is homogeneous and isotropic, the 
coordinates can be chosen in such a way that the base 
plane for which X 2 = 0 passes through Xa^b- Then for 
given Xa^Bi the correlator Vab is determined from the 
1+1-dimensional capillary wave held theory coupled to 
the 2+1-dimensional fermions, 

( 23 ) 

with the action S = SeF Sew + ‘^int- 

We emphasize here the subtlety of dimensionality in 
our approach: (1) We are concerned with the long 
time dynamics of the system, and the system is ho¬ 
mogeneous and isotropic on average at this time scale. 
Hence the correlator Vab is defined in 2+1-dimensions, 
and depends only on the separation between Xa^b^ i-e. 
Vab = V{\Xa — Xb\)‘ We can then Fourier transform 
Vab from the 2+1-dimensional space-time to the mo¬ 
mentum/frequency domain. (2) To calculate Vab^ we 
need to specify the two points Xa,b^ and sum over all 
interfaces that pass through them. We will use the cap¬ 
illary wave formalism for this step, that is, we sum over 
low energy fluctuations around the base plane that passes 
through Xa^B‘ Such capillary wave fluctuations are en¬ 
coded in a 1+1-dimensional held theory, in terms of the 
transverse mode of the string, i.e. the height function 
h{t,x). (3) Fermions couple to the interfaces and af¬ 

fect their dynamics, causing Landau damping. Since 
fermions live in 2+1-dimensional space-time, we will need 
to project their effect to the 1+1-dimensional interface. 

We want to compute the Fourier transform of Vab 
to the momentum/frequency domain, which can be ob¬ 
tained from the correlator P(r, cj) 

V(q,u) = j d?re-^‘^-^V(r,w). (24) 

The correlator V{r^uj)^ where r = 1^^ — is 

related to the capillary wave correlator Dh{k^uj) = 
{h{k^uj)h{—k, —uj)) via 

V{r,uj)= [ dkF^^E^k^Dh{k,uj). (25) 
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After integrating over the angular coordinate, one ob- The fermion self-energy then reads 
tains the following relation between V{q^ uj) and the cap¬ 
illary wave correlator ^ iArcSinh(cj/'z;sg'o)- (32) 


V{q,uj) 


j drrJo{qr) j i 


dke^^^E^k^Dh{k,uj), (26) 


with the Bessel function J. 

Once V{q,uj) is known, the one-loop fermion self¬ 
energy can be obtained from 


E{k, uj) = J d?qdrtV{q, Q)G{k + gr, cj + f 2 ). (27) 


Since the bosonic mode is much slower than the fermions, 
the momentum integration can be factorized into two 
parts f dq± f dgy, with q± representing the momentum 
component transverse to the Fermi surface, i.e. the fast 
fermion modes, gy the parallel component, i.e. the slow 
modes. Restricting k to be near the Fermi surface, and 
carrying out the q± integral, one obtains |47l48j 

E(u;) ^ iA^ f dftd{ft)sign{uj ^ ft). (28) 


where the bosonic correlator appears through the 
momentum-integrated “local” form 


rqo 

d{uj)= dqiiV{qii,qj_ =0,u)), (29) 

Jo 

with the momentum cutoff q^ ^ kp. 

We have presented above a general procedure to cal¬ 
culate the fermion self-energy from the 1 + 1 -d capillary 
wave correlator Dh{k,uj). We now proceed to consider 
different models of the interface fluctuations that give 
specific forms of Dh{k,uj). 


A. Fermion self-energy for perturbative models 


After analytic continuation to real frequency, one can 
see that the imaginary part of the retarded self-energy 
ImSi^(co’) = 0 for cj < Vsqo. (2) When the dynamics is 
controlled by Landau damping, i.e. aq^ + 7 |cc;|g^, 

one has 


V{q, w) = d{uj) -TTTj- 

The self-energy then reads 

Il(cj) ^ isgn(cj)log . (34) 

Analytic continuation to real frequency yields 
lmER{uj) log At low frequencies 

UJ ^ cr/ 7 , ImSi^(u;) ^ which is of the Fermi liquid 
form. Hence for large interface tension, the fermions 
remain in the Fermi liquid phase. 

When the interface tension vanishes, i.e. a ^ 0, the 
interface fluctuates much more violently. In this case, 
the dynamics is controlled by the Landau damping term, 
with 


Dh{q,oj) 


1 

Kq'^ + j\u:\q^' 


The resulting string correlator is 


V{q,uj) r-. 


1 

{Kq^ + 


(35) 


(36) 


The dynamical exponent which relates the temporal di¬ 
rection to the spatial direction in cj ^ is z = 2. The 
momentum-integrated local string correlator reads 


One type of models that take a more perturbative point 
of view is to start from Eq.(13), where one considers the 
interface to have a tension a, a curvature stiffness k, and 
a velocity Vs relating frequency and momentum. Cou¬ 
pling to fermions induces extra dynamics for the interface 
fluctuations. Including such Landau damping effects, the 
capillary wave correlator reads (see Appendix for the cal¬ 
culation) 


Dh{q,oj) 


1 

cr(g2 _|_ Vs^Lu'^) + K{q^ + Va‘^u)‘^Y + 'y\uj\q^ ’ 

(30) 


with a new parameter 7 oc A^. 

When the interface tension is large, the curvature term 
can be neglected. We consider two limiting cases. (1) 
When the Landau damping term is small, the dynamics 
is relativistic, i.e. 'Zi cF{q^ -\-vj‘^uj‘^). The local string 
correlator is of the form 


d{uj) ^ 


Vsqo 
v/w 2 + 


( 31 ) 


d{i^) ^ A, (37) 

|w| 

and the fermion self-energy is 

Il(u;) ^ isgn(cj) log |u;|, (38) 

which diverges as cj ^ 0 . 

We can compare our result with other types of fluctu¬ 
ations. For example, quantum critical fluctuations as¬ 
sociated with 2d antiferromagnetic QCPs give rise to 
E{uj) ^ ^JuJ [49ll5T] . which is much milder than our re¬ 
sult of a divergent self-energy. Such ^JuJ form of the self¬ 
energy was also obtained in the random infinite-range 
interacting model [52ll54] . Similarly, quantum critical 
fluctuations associated with 2 d ferromagnetic QCPs lead 
to S(co’) ^ uj^l‘^ [55]. It is surprising that the innocent¬ 
looking capillary wave field theory should have such a 
dramatic effect on the fermions. In fact, the capillary 
wave field theory can be expressed in a more familiar 
form by defining 0 = dh/dx. Then the tension term is 
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just the mass term of the bosonic field and the cur¬ 
vature term the gradient term Setting a ^ 0, 

one obtains a massless bosonic field coupled to fermions 
that is well studied and no such singular effect was found. 
The enhancement of the self-energy correction in our case 
comes from the fact that the interface fluctuations reside 
in one dimension lower than that of the fermions. It is 
well known that fluctuations are more severe in lower 
dimensions. 


B. Fermion self-energy for Luttinger-liquid type 
models 


To put the above results in a broader perspective, let us 
consider another type of model that takes a more scaling 
point of view. From the action as shown in Eq.(14), one 
obtains the capillary wave correlator 


Dh{q,i^) 


(g2 + t,7^a;2)i+C 


(39) 


The resulting local string propagator is 

and the fermion self-energy is of the form 

S(u;) ^ isgn(ci;)|ci;|^“^^. (41) 


The self-energy diverges for rough surfaces with > 1/2. 
For ( = 1/2, where d{uj) ^ 1/|^|, the self-energy is 
F(cc;) ^ isgn(co’) log |cc;|, which diverges logarithmically. 
This case corresponds to the above perturbative model 
with vanishing tension a = 0. 

We can take an even more coarse-grained scaling point 
of view to regard the local string propagator d{uj) as 
having a certain scaling form, and use that as input 
for the calculation of the fermion self-energy. It is in¬ 
teresting to note that the marginal fermi liquid with 
F(u;) ^ iuj\og{ujo/uj) corresponds to d{uj) r\j log|w|, i.e. 

C^o. 


Due to the singular frequency dependence, the momen¬ 
tum dependence becomes negligible, and the full fermion 
Green’s function becomes local. Furthermore, the origi¬ 
nal pole of free Fermi gas Green’s function at cj = 0 is 
converted to a zero by the divergent self-energy, which 
indicates that propagating single electron states cease 
to exist. The Landau Fermi liquid picture completely 
breaks down, and we obtain a non-Fermi liquid, i.e. 
strange metal state. The absence of propagating sin¬ 
gle particle states in our system is intuitively clear: at 
the critical point, where the interface tension vanishes, 
the interface forms sponge-like structures at large scales 
(see Figj^d)). Hence the electrons scatter so much with 
the interfaces on their way of propagation that they can 
no longer propagate at large spatial/temporal scales, and 
are effective confined. 

It is interesting to note that, the conversion of poles of 
the fermion propagator to zeros also occurs in the 1 + 1-d 
large-gauge theory, which is a toy model for quark con¬ 
finement [56]. There such behavior is regarded as a clear 
signature of the absence of physical single quark states. 
The significance of the zeros of the Green’s function and 
in particular their relation with the Luttinger theorem 
have been discussed in [57158] . 

With the absence of propagating single particle states 
at low energies, the resistivity of such non-Fermi liquid 
metals is no longer controlled by the quasiparticle decay 
rate as in a Fermi liquid metal, but by the decay rate 
of the total momentum. Since the system equilibrates 
very rapidly, at scales shorter than that of the breaking 
of momentum conservation due to, for example, impurity 
and Umklapp scattering, the transport properties of the 
system can be described by hydrodynamics [59160] . The 
temperature dependence of resistivity is then determined 
by the scaling dimension of energy density and charge 
density operators noil]. One way to compute the resis¬ 
tivity for systems near the MIT is to employ the random 
resistor network |32| . The quantum-critical sponge phase 
identified here corresponds to the spin liquid phase of a 
frustrated Ising model. It would be interesting to cal¬ 
culate the resistivity of the resistor network defined on 
a geometrically frustrated lattice. We leave this for a 
future investigation. 


C. Zero of the Green’s function 


We found above that the interface fluctuations can 
give rise to a divergent fermion self-energy. The fermion 
Green’s function G{k^uj) = is then domi¬ 

nated by the self-energy term at low frequencies. For 
the perturbative model with a = 0, one has 


V. INTERFACE FLUCTUATION MEDIATED 
PAIRING 


The interface fluctuations induce a density-density in¬ 
teraction among the fermions as shown in Eq.(20). It can 
be written more explicitly as 


GiKu) ~ ( 42 ) 

logical 

For the Luttinger-liquid type model with > 1/2, one 
has 

G{k,uj) ^ isgn(cj)|u;|^^“^. (43) 


^ J dti /-^■/ dh j d'^r 2 !){r,t) 

^Uri,ti)^pa{ri,ti)ipl{r2,t2)^ppir2,t2) (44) 


with the bosonic string correlator P, and r = ri — r 2 , t = 
— ^ 2 - Here we have restored the spin indices a, [3 =t, 














and the summation over them is kept implicit to simplify 
the notation. 

As is familiar from the case of phonon mediated pair¬ 
ing, such density-density interaction can give rise to 
Cooper pairing in the spin singlet channel. We note 
the difference with models employing homogeneous ap¬ 
proaches |21l62j . where a repulsive current-current in¬ 
teraction is obtained between two fermions with oppo¬ 
site momentum that pair [63]. In our case, the local 
electric field resulting from inversion symmetry breaking 
around the interface induces a direct coupling between 
the fermion density and the bosonic mode, which leads 
to an attractive density-density interaction among the 
fermions. However, as shown above, the interface fluctu¬ 
ations also generate self-energy corrections for fermions, 
which lead to incoherence, and are actually detrimental 
for superconductivity. Indeed for a = 0 in the pertur¬ 
bative models, and for ( > 1/2 for the Luttinger-liquid 
type models, the self-energy diverges. The pole in the sin¬ 
gle electron Green’s function is converted to a zero, and 
there is no physical single electron state. However, the 
existence of Cooper pairs is associated with the pole of 
the two electron Green’s function in the particle-particle 
channel. To see whether such a pole exists at the one- 
loop level, we can use the Eliashberg formalism, which 
takes into account the competing effects of self-energy 
and density-density interaction. The Eliashberg formal¬ 
ism has been applied to systems of fermions coupled to 
quantum critical fluctuations in j47l48l51l64l65j . 

The density-density interaction can be decomposed in 
the pairing channel. Keeping terms relevant for zero 
total-momentum spin-singlet pairing, one has 


- -A^ J <fk J dco J J du'V{k -k',uj- u') 

x\k,u})x{k',uj'), (45) 


with the spin-singlet pairing operator = 

-^{'>PkLo-fi’-k-Lo,i - i’kLoii’-k-Lo,-f)- The superconduct- 
ing gap can then be defined from 


A(fc,w„) = A^T ^ -D{k - -0(x(fe',w;,))(.46) 

1^'m 


Since the interface fluctuation mediated interaction is 
isotropic, pairing is in the 5-wave channel. Other de¬ 
grees of freedom, e.g. spin fluctuations [66], need to be 
invoked to generate unconventional, e.g. d-wave, pairing. 

With the fermions moving much faster than the 
bosons, momentum can integrated out in the Eliash¬ 
berg formalism, leaving only the frequency dependence 
explicit. Near the superconducting transition, one has 
the coupled equations for the self-energy (EqlH 

41) and the pairing vertex (see. e.g. [48]) 


= ttt y] 


d{uJjn ^n) 

\^m \ 




(47) 


with the pairing vertex 4>(ci;^) = A(u;^)Z(co’^), and the 
renormalization factor Z(uJn) = 1 + Pairing 


Tc/Qo 



FIG. 4: (Color online) The pairing temperature as function 
of the critical exponent ( in the bosonic propagator. 


occurs when the gap equation (Eq.(47)) have a solution. 
In particular, for a quantum critical boson with the local 
correlator d(u;) = (flo/l^l)^^, where 0 < C < Eq.(47) 

always has a solution. Since flo is the only dimensionful 
parameter here, the resulting pairing temperature Tc is 
set by flo, with |4VI48l5H 


Te = llo/(C)- (48) 

The coefBcient f(() decreases with increasing ( (see 
FigQ, implying a stronger effect of incoherence. 

The above calculations imply a novel picture for 
strange metals: while there are no propagating single par¬ 
ticle excitations, there can be Gooper pairs. Intuitively 
one expects fermions confined in the sponge structures to 
be unstable, since their zero point energy is significantly 
boosted by the confinement. Pairing is a plausible way to 
release such energy. In fact, superconductivity has been 
widely observed near Eermi-surface-changing MIT in ma¬ 
terials like cuprates mm, organic conductors [22|23], 
iron pnictides |24l25j and heavy fermions |26l27j . How¬ 
ever one needs to be cautious with the above perturba¬ 
tive treatment of pairing. We discuss below a more local 
picture of pairing in the fluctuating sponge structure. 

At time scales where the interface can be treated as 
quasistatic, the electrons are localized to cells of size 
due to the inhomogeneous background. Then Gooper 
pairing occurs locally. In 2d the binding energy of two 
electrons scales with the size of the system as ^ 
1/(1/^ log!/) (see Appendix H). With L ^ pairing 
is largely enhanced by confining electrons in the sponge 
structures. This effect is in some sense the real space 
counterpart of Gooper’s mechanism of pairing [67]. The 
celebrated Gooper instability relies on the existence of the 
Eermi surface, and electron binding is enhanced due to 
Pauli blocking in momentum space. In our case, electron 
binding is enhanced due to the temporary confinement of 
electrons in real space. The occurence of superconduc¬ 
tivity also requires phase coherence. In the fluctuating 
sponge strucutre, electrons can first pair locally, taking 
advantage of the temporary confinement of the local ge- 
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s 



FIG. 5: (Color online) The entropy as function of temperature 
for three different characteristic interaction scales. 


A ^ 0.73). Since is the only dimensionful param¬ 
eter here (besides the cutoff), it is not surprising that 
T* is on the order of The negative entropy implies 
that the strange metal phase becomes thermodynami¬ 
cally unstable for temperatures T < T* ^ Indeed, 
as shown above, pairing sets in at the same temperature 
scale Tc ^ flo- The concurrence of the two temperature 
scales T* ^ Tc originates from the fact both the self¬ 
energy and the pairing interaction have the same char¬ 
acteristic energy scale Qq ^ f P, which is determined 
by the coupling strength A and the bosonic propagator 

V. The resulting picture is that the strange metal phase 
is intrinsically unstable, and it is forced thermodynami¬ 
cally to be replaced by a symmery-broken phase at low 
temperatures. 


ometry. And then, at larger time scales, where the ge¬ 
ometric fluctuations come into action, the whole system 
becomes phase coherent. Hence pairing and phase co¬ 
herence have different time scales, and correspondingly 
different energy scales. A more in-depth investigation of 
this pairing mechanism is left for a future publication. 


VI. THERMODYNAMIC INSTABILITY 


We have shown above that the interface fluctuations 
induce pairing instability among the electrons. Here we 
will show that there is indeed a deeper reason for the 
occurrence of pairing: the strange metal phase is ther¬ 
modynamically unstable below certain temperature T*, 
and hence some new phase must set in at such a tem¬ 
perature. The corresponding temperature scale is on the 
order of the characteristic interaction energy scale flo, 
which, as shown above, also determines the pairing scale, 
i.e. T* ^ Tc ^ ^ 0 * 

The entropy of the system can be determined from the 
relation [68] 


S 

V 


-1 

iriT 



[logG'^,c^ -logG^,^](,49) 


where V is the volume of the system, and / is the Fermi 
function. and ^ are the retarded and advanced 

fermion Green’s functions, and are obtained from the 
Matsubara Green’s function through analytic continua¬ 
tion. For the self-energy of the form as shown in Eq.(38), 
they read 


qR,a ^ _ 1 _ 

- f^^osgn(w) zbifiolog]^’ 


(50) 


with the characteristic interaction scale Qq, and the cut¬ 
off A. 

The results are shown in Fig. We find that the 
entropy decreases with lowering temperature, and it be¬ 
comes negative at T* = with A of order unity (here 


VII. CONCLUSIONS 

We have proposed a new framework for non-Fermi liq¬ 
uid metals near a Mott metal-insulator-transition with 
fermions operating in sponge-like structures formed by 
strongly fluctuating interfaces. The geometric fluctu¬ 
ations of the interface give rise to anomalous scalings 
in the fermion self-energy that are qualitatively differ¬ 
ent from those of Landau Fermi liquids. Furthermore 
the geometric fluctuations induce pairing instability in 
the absence of external bosonic glues. The resulting pic¬ 
ture of a strange metal phase at high temperatures and 
a low temperature superconducting phase is in qualita¬ 
tive agreement with the experimental phase diagram of 
cuprates, iron pnictides, organic conductors and heavy 
fermions. 

By directly incorporating quantum geometric fluctua¬ 
tions into our framework, our result constitutes a qualita¬ 
tively new step going beyond the homogeneous paradigm 
of many body systems. The fluctuating geometry is not 
a perturbation to the system, but instead provides the 
template for all other degrees of freedom. It will be in¬ 
teresting to study other processes, e.g. spin-fluctuation 
mediated pairing [69|, in such a fluctuating background. 

A crucial question is how to distinguish between the ef¬ 
fects of geometric fluctuations and fluctuations of other 
degrees of freedom. One strategy is to look for strange 
metal behavior away from the symmetry-breaking QCPs. 
In fact, strange metal behavior was observed in Ir- and 
Ge-substituted YbRh 2 Si 2 [THTfl] and ;5-YbAlB4 under 
pressure [72] in regions associated with a Mott transi¬ 
tion away from the magnetic QCPs. Quantum geometric 
fluctuations may play important roles in these systems. 

One further extension of our work is to consider fluctu¬ 
ating interfaces involving topologically nontrivial phases, 
for which there can be gapless fermion modes emerging 
at the interface. A concrete example is a heterostructure 
of a 3d topological insulator with a magnetic insulator 
near the ferromagnetic transition. When approaching 
the transtion from the ordered side, magnetic domains 
will be created. A gapless chiral fermion mode is then 













10 


induced on the topological insulator at the location of 
the domain wall The domain walls are then het¬ 

erotic strings m in the sense that the fermionic mode 
is chiral while the bosonic modes have both left- and the 
right-moving components. 
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Appendix I: Landau damping 


Keeping only the leading power law dependence, the re¬ 
sult is of the form 


Uhik,Lo) ='y\Lo\k'^, (58) 

with the prefactor 7 cx X'^E^. We need to keep in mind 
that 7 depends on the IR cutoff. The capillary wave 
correlator thus reads 


Dh{k,uj) 


1 

(59) 


We also note that with the full Green’s function 
G{k^uj) ^ isgn(u;)/log |cj|, the fermion bubble gives 



sgn(cj) sgn(cj + fl) 
^ log |cj| log Icj + Q\ 


- \n\. 


(60) 


The leading power dependence of Ilh{k,uj) is then the 
same as using the free fermion Green’s function, and the 
above form of Dh(k^uj) is self-consistent. 


Goupling to fermions leads to Landau damping in the 
capillary wave correlator. One has the Dyson equation 

= (51) 

with the free part 

[d7]-i = a{e + + K{e + (52) 

and the self energy correction 11^ coming from coupling 
to fermions. The 2+1-dimensional fermion bubble 

Ilo{q,Q.) = j -^^^Go{k,u})Go{k + q,u} + Q.) (53) 

is of the form 


no(q,fl)-|fl|/</. (54) 

Eourier transforming to 2+1-dimensional real space 
yields 

no(r, n)= [ ^e^i-Uoiq, Q) ~ (55) 

J (27r)^ r 

Since the system is homogeneous and isotropic on aver¬ 
age, no(r, 0 ) = no(r, O) determines the self energy in 
the 1+1-dimensional field theory. With the fermion den¬ 
sity coupling directly to one obtains the self energy 
correction for {dhdh) as 

Uoh{k,oj) = X^E^ [ dre-''=’'no(r,w) - -( 70 +log |fc|)|!d|, 

(56) 

with Euler’s constant 70 . We note that there is a loga¬ 
rithmic divergence at small momenta, and an IR cutoff 
needs to be imposed. The self energy correction for the 
{hh) correlator is then 


(57) 


Appendix II: Cooper pairing in a confining geometry 


Consider two particles in a harmonic trap interacting 
via a pointlike force with Hamiltonian H = Hq + iLi, 
where [76] 

Hi = (62) 

with the scattering length uq, and the regularized 6 func¬ 
tion potential Jreg(^)- Such a point-like potential is a 
good approximation for uq ^ T, where L = y^H/(mujo) 
is the characteristic length scale of the harmonic trap. 
The characteristic energy scale of the harmonic trap is 
El = kwo = h^t{mL^). 

In 3d, the binding energy is related to the size of the 
trap via 



When the trap size is on the order of the scattering 
length, L ^ ao, the binding energy is ^ h‘^/{maQ). 
El) increases rapidly as L decreases, and decays rapidly 
as L increases. Mathematically bound states always exist 
in the presence of a trap. But practically the binding en¬ 
ergy is too small to lead to any observal effect for L ^ ao. 
At large L, one has 


E, 


k?aQ 

mL^ 


In 2d, the binding energy is given by m 


(64) 


ip 


Eb 

2El 


= log 



Jlh{k,Lo) = fc^na/i(fc,w). 


(65) 











11 


where is the digamma function. For x ^ 0^ pj{x) 
—7 — Ijz. At large L, one has 


Eb 


mL‘^ log(I/^/aQ) 


( 66 ) 
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